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We investigate the approach to the quantum critical point of a Pomeranchuk 
instability from the symmetric, disordered side of the phase diagram. In the 
low-temperature limit, a Fermi liquid description of the metal is possible and 
becomes exact for T — > 0. We discuss in detail which features of the approach 
to quantum criticality can be captured within Fermi liquid theory and which 
are outside of its scope. 
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1. Introduction 

Fermi liquid theory is the basis of our understanding of metals. In 
recent years, quantum phase transitions in metals have been widely studied 
mainly motivated by the experimentally observed (apparent) violation of 
Fermi liquid behavior for a wide variety of system^^J^. While a Fermi 
liquid description may break down directly at a quantum critical point, it 
is expected to be fully valid arbitrarily close to this singular point in the 
low-temperature limit, T ^ (at least in the symmetric phase). Therefore 
we want to explore in this paper the Fermi liquid description of this regime. 

From the point of view of Fermi liquid theory, the simplest and most 
straightforward instability is one where the Fermi surface is deformed but 
does not change its topology. Such PomeranchulP instabilities have been 
widely studied in the last few years I ^ I ? I 7 I ?I9 I 19I1J^ | 12 |_ Experimentally, spon- 
taneous deformations of correlated quantum systems, e.g. of heavy nuclei 
or of crystalline lattices, have been well-known for decades. More recently, 
purely electronic nematic transitions have been suggested to be relevant for 
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systems ranging from high-temperature superconductors and 2-d Hubbard 
models l^'^^^^, Sr3Ru207 close to its metamagnetic transitiori^to clean two- 
dimensional electron gases in high Landau levelPl. Nevertheless, we are not 
aware of experiments on quantum critical properties of a nematic transition 
e.g. from a cubic to a tetragonal phase of a metal. An exciting option 
for the (near) future is the controlled realization of a nematic transition of 
Fermions in cold atom systems using Feshbach resonances in channels with 
higher angular momentum. In these systems it should be possible to excite 
and measure directly the relevant collective excitations by studying defor- 
mations of the atom cloud in response to changes in the trapping potential. 

Probably the most intriguing aspects of nematic phases in isotropic 
systems is that deep in the ordered phase Fermi liquid theory is not valid as 
pointed out by Oganesyan, Kivelson and FradkiiP. As explained above, we 
will, however, focus our attention on the low-temperature behavior of the 
phase where the symmetry is not broken. We first briefly outline the Fermi 
liquid setup, then discuss the collective modes, the self-energy and finally in 
the concluding section how Fermi liquid renormalizations couple back to the 
order parameter dynamics. 



2. Fermi liquid model 

In the Fermi liquid phase of a metal fermionic quasiparticles are the 
dominant low-energy excitations. At low temperatures or, more generally, in 
weakly excited states the number of quasiparticles or holes is small, providing 
a small expansion parameter. At the same time the quasiparticle interaction 
may be strong, as measured in terms of dimensionless Landau parameters 
. For an isotropic system, which we consider in the following, the single 
quasiparticle energy 

efc = VF{k - kp) (1) 

where kp is the Fermi momentum and vp = kp/rn* the Fermi velocity, is 
isotropic, and the Fermi liquid interaction -Ffco-fc'o-' may be expandecf^ in 
terms of spherical harmonics, 

1 °° 

fkak'.' = ^Yl m-k'm + Flaa'] (2) 

^ i=o 

where a = ±1 are the spin quantum numbers, Np the quasiparticle density 
of states and Pe{x) are in d = 3 dimensions the Legendre polynomials, 
Pi{k ■ k') = {A'K/{2l+l)}T:myim(^)yim{ki), while in d = 2 dimensions Pt{k ■ 
k') = cos£{(j) — (j)') (with (/), (p' the angles of k, k' with the x-axis). A change in 
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the quasiparticle occupation number Jn^cr leads to the quasiparticle energy 
change 

S^ka = y ^ fkak'a'^nk'a' (3) 
fc'cr' 

where V is the volume of the system. A spatially non-uniform change of 
occupation number with Fourier component q, 6nka{Q), will lead to a corre- 
sponding energy change 5eka{Q)- For sufficiently slow spatial variation such 
that q <^ kp, the Fermi wave number, the dependence of fkak'a' on q may 
be neglected (it gives, however, a contribution^ to the q dependence of the 
order-parameter susceptiblity, discussed below). 

On the microscopic level we will adopt the model Hamiltoniar^ 

I ' 

H = '^eknka + fkcTk'a'nka{q)nk'a'{-q) , (4) 

ka k,k',q 

where the prime on the summation sign indicates that the g-summation is 
restricted to q < qo <^ kp- 

In a crystal lattice the full rotation symmetry is reduced to the point 
group symmetry. This implies that the Landau parameters do not only 
depend on i, but also on m. 



3. Fluctuation propagator and Pomeranchuk instability 

The response of the system to an external field uj) = J2em ^em{k) 

[^^Im'^ + '^^^em'^]^ where Yim{k) are the spherical harmonics in d = 3 (in 
d = 2 dimensions the expansion is in terms of sml(j),cos i(p), is expressed 
for each channel {i, m,J = s, a) in terms of the generalized susceptibility 

xUq,-) = ^. (5) 

We assume the vector q to be oriented along the z-axis in d=3 and along 
the X-axis in d=2. Within RPA xim given by 

" l + F/Ue,m/2NF{2e + l) 
where Hi^m is the susceptibility in the absence of Fermi liquid interactions. 



UJ + iO- €k+g/2 + ek-q/2 
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Expanding this expression in small q and {lo/vfq) one finds at low temper- 
atures 

Ue,m{q,uj) = 2Np[l-eo^^,^q^+P'emi^,q)+il'imi^,q)+0{q^)] (8) 

where Np is the density of states (per spin) of quasiparticles at the Fermi 
level. The functions and 7^^, discussed below, depend on whether 
the considered mode is "even" or "odd", which means whether in d = 3 
dimensions ^ + m is even or odd and in d = 2 whether the mode is cos(2n(/)) 
or sin({2n + l}^) or whether it is sin(2n0) or cos({2n + !}</>). 

At this point one has to investigate which properties of the dynamic 
susceptibility X'emil^'^) calculated within Fermi liquid theory. Both 

the to = response and the damping of the collective modes in the limit 
a; — > are low-energy properties which can be described exactly using Fermi 
liquid theory. However, this is not the case for the momentum dependence 
which gets contributions from high-energies, i.e. from the incoherent part 
of the excitation spectrum. Especially, close to the critical point when the 
quasiparticle weight vanishes (see below), one can expect that the response is 
completely dominated by the incoherent part and a Fermi liquid calculation 
is impossible (further contributions also arise from momentum dependencies 
of the screened interactions). However, it is simple to give an order-of- 
magnitude estimate of this term by assuming that only renormalizations 
of order 1 occur for this high-energy quantity. Taking into account that 
Np (X 1/Z , where Z = m/m* is the quasiparticle weight factor, and a 
relative factor of between the polarization 11 of quasiparticles Q and 
electrons, we estimate 

fully consistent with results in the literatur d^ l ^ l ^ l 

In contrast, the damping term, , can be calculated rigorously from 
Fermi liquid theory, yielding in d = 3 



7^^(w, q) = Att Yi^icos 9 = , </.) (10) 



vpq 



vpq 



The coefficient of in this expression is finite in the limit cj/vpq — > 0, or of 
order (uj/vpq)'^, depending on whether £+m is even or odd. This dependence 
arises because only quasiparticles with velocity perpendicular to q, i.e. with 
angle 9 ~ 7r/2, contribute to Landau damping due to momentum and energy 
conservation. This has the consequence that even and odd modes have 
very different dynamical exponents. As for the even modes the Fermi liquid 
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renormalizations cancel between momentum and frequency dependence, one 
can directly read off the dynamical critical exponent ^even ~ 

3. To discuss 

Zodd we first have to calculate the scaling of m* /m upon approaching the 
QCP which is done in the next section. 

The real part and the imaginary part of Hi^miQ,^) are connected by a 
Kramers-Kronig relation. In case of the even modes the imaginary part does 
not give rise to a singular contribution to the real part, i.e. fi'lm"" ™a,y be 
neglected. However, for the odd modes, one finds a contribution to the real 
part 

2 



a/odd 



{uj,q) = 8iT Yijn{cos6 



CO 



vpq 



\VFqJ 



(11) 



In d = 2 dimensions the damping term is similarly suppressed whenever 
one of the functions sin ^0, cos vanishes at (j) = ti/2. Note that the suppres- 
sion of damping is only effective for an isotropic system. In a crystal lattice 
the eigenfunctions are pinned to the lattice directions, and the condition of 
reduced damping can be only satisfied for special directions of q. 

In the limit of small q and ujjvpq we get finally 



2NF{m/m*) 



(12) 



where the correlation length the damping term 
term are defined as 



and the inertial 



7emi^,q) = 

Pemi^^q) = 
in d = 3 dimensions and 



1 + 



2^ + 1 



m 



(13) 



m 
m* 
m 
m* 



a 



io/vpq 



(^imi^/^Fqf 



\2 



,i + m even 
,i + m odd 



: + m even 
bj^{uj/vFqf ,i + m odd 



^2 

iimi^^q) 
Pimi^^q) 



m* 



(14) 



m 
m* 
m 
m* 



a 



im '^hpq 
(^Imi^/vpqf 


Kmi^/vFq? 



, for cos £(l) 
, for sin £(j) 

, for cos £cj) . 
, for sin £6 , 



£ even; sin £4) , £ odd 
? even; cos £(p , £ odd 

' even; sin^c/) , £ odd 
' even; cos £(h , £ odd 



in d = 2 dimensions with numerical constants a^^' ' . As we will show below, 
for even £ + m our results for x are fully consistent with the literatur e! '"^ I ^ l ^ l 
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Fermi liquid renormalization effects are much more important for the odd 
modes (see below), which have either not been investigated^ (they are only 
relevant for uniform systems) or even when considerecP' the Fermi liquid 
corrections have not been studied. 

The uniform static susceptibility x/m.(0;0) is seen to diverge when 

F/ ^ -{2£ + l), d=3 (15) 
F/ ^ -1, d=2 (16) 

signalling the Pomeranchuk instabilitjEl. For £ = 0, J = a this is nothing but 
the familiar ferromagnetic instability. We will, however, only consider non- 
magnetic instabilities as ferromagnetic transitions, where the Fermi surface 
splits, turn out to be substantially more singular (Ref. 8 and references 
therein). 

Near the quantum critical point the fluctuations induce an effective 
dynamical quasiparticle interaction of the form 




whose effects are studied in the next section. 

Note that for d = 3 and even £ the transition is generically of first order 
due to the presence of cubic terms in the Ginzburg Landau expansion which 
are absent for odd £ in d = 3 and for arbitrary £ in d = 2. Also a transition 
in the charge channel with £ = 1 is not possible for a Galilean invariant 
systems as the ordered state with 1 = 1 describes a Galilei transformation 
which always costs energjW]. In the following we will tacitly assume that a 
first order transition does not take place (or is sufficiently weak). 



4. Self-energy 



The contribution to the quasiparticle self-energy induced by critical fluc- 
tuations in channel £ takes the form 



Im S(k, w) 



' /o J Im{xL(g,u;'-iO)}5(a;' + u;-ek+q)(18) 



n 



where nP, and fP are Bose and Fermi distribution functions. We recall 
that the spherical harmonics are defined with respect to the axis q. Since 
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Im Ximil^^) is isotropic in q, the angular q-integration removes the depen- 
dence on k , such that S(k, w) is a function of |k| only, as expected for 
an isotropic system. At low temperatures, and small uj , k the dominant 
contribution to the integral comes from small frequencies and by the energy 
conserving delta function, small k • ^ . As discussed above, the quantity 

2 



47r 



Yf 



in this limit will be a constant for even modes and equal to 

x^m(k • q)'^ for odd modes. 

The uj'- integration in Eq. (jlHj) may be readily performed. As for the 
q-integration, we use the decomposition q = g^k -|- qt , where qt is a (d — 1)- 
dimensional vector in the plane perpendicular to k and k = k/|k|. The 
quasiparticle energy ek+q may be expanded in qr as 



ek+q = {k + qr)vF 



(19) 



where vp is the Fermi velocity. Then, 



Im S(k, io) 



2Nf 



dqr f d'^qt 
27r J (27r)2 



n^{ek+q - + /°(efc+g) Im{xL{q , 



At zero temperature the distribution functions reduce to step functions, e.g. 
for w > , 



-1 



0<e<a;or -u < e <0 



, else 

and the contribution of even modes im to S is then 



Im uj) 



{F^f ^ dqr f _At VimlioJ, q) 



2vr J {'i^fiilC + iW-Pi^.qW 



(21) 



r 



(22) 

where u = uj — vpik + qr) . Here we dropped the indices £m of 7^^, and 
^Q£rn introduced the correlation length ^ of mode im by 



,m 



^ =^o,imi—) l + Fe'/{2i + l) 



-1 



(23) 



It is seen that the main contribution to the q-integral comes from regions 
where qt ^ qr, and therefore we may approximate g ~ g*. 

If ^ -|- m is even, and 7 = auj/vpq, /3 = 0, one finds with the new 
dimensionless variables ijr = vp^qr + k)/uj and qt = qtjquj, where q^ = 
{auj /v^p^'^y/'^ , and v^p is the bare Fermi velocity. 
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2Nf a^^ Jo 27r J (27r)'i-i ^2(^2 ^ -sp + _ ^-^p 

(24) 

where ( = {quj^,) ^. Here we have generalized the expression to arbitrary 
dimension d. 

Now the integration on ijr may be performed to give 

where Q = {q'f + C^)qt. 

For small w, <C = vp^'o/i^a., we have C 1 Q ~ C'^Qt- Using 
that Q ^ 1 in most of the integration regime, we finally get 

Similarly one finds for a two-dimensional system 



ImS(^,.) = ^Al|^(^)- Ini^) e (27) 



consistent with results of Dell'Anna and Metznef ^ if one takes into account 
that in H27() we have calculated the self-energy of quasiparticles rather than 
the one of electrons which is a factor \ — larger. 

Within Fermi liquid theory one can calculate exactly the quasi-particle 
decay rates in the limit of vanishing quasi-particle energy as higher-order 
contributions beyond Eq. (PU]) are supressed by phase-space factors. How- 
ever, it is generally not possible to determine effective masses and the quasi- 
particle weight Z which are an input to rather than an output of Fermi 
liquid theory. Close to the quantum critical point, however, the singular 
corrections to real- and imaginary part of the self energy come from the 
same critical fluctuations. Therefore we can proceed by obtaining the real 
part of S from Eq. by replacing £n(l + ^) = Im{+i£n(l + ^)} by 

Re{i^n(l + ^)} = — arctan ^. It follows that 

ReS(fe,u;) = -^^^ln(C/U, ^ = 3 (28) 

which is still a quantity obtained from a quasiparticle description. To cal- 
culate Z ~ m/m* we need the self-energy of the original electrons Tf^\k,ijj) 
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(rather than of the quasiparticles). But this quantity can be obtained using 
S^'(A;,a;) = (m*/m)S(/c, w), to give 

m* 5i?eE^' _ 1 / 1 \ 2 , 



Only a logarithmic enhancement arises for ^ — > oo in d = 3. 
In two dimensions, however, we find 

ReE(A:,u;) = --^ i (30) 

^ ' ' 327r INp v%il ^ ' 



and hence 



^ = i + ^(feFCor'(^/)'(^/u (31) 



which is linearly diverging as ^ — > oo. 

In the case of odd £ + m, where the additional contribution /3 appears 
we may consider the limit 7 ^ 0, as 7 is a subleading term compared to /3, 
so that Imx reduces to a delta function. One finds 

ImE(fc,a;) = ^{^) 2^ 7 (^^^^^i^ ^^'^ 

^,2,2 _^( A.).] (33) 

At the Fermi surface (A: = 0), performing the qr - integration with the help 
of the delta function and defining the dimensionless variable p = qt^ one 
finds 

""^(*'") = 4(l)^ ^--^f ^ C (34) 

where = {m* /m){u) /v^p^^f'Sf^ , and Via is the surface area of the unit 
sphere in d dimensions. We have used that in the limit of small frequencies 
we may neglect compared to 1. In d = 3 dimensions we find 

In two dimensions one observes that the p integral is cutoff at the lower 
end by the requirement qt > qr , yielding p > Pc = (^)i^)- For odd i + m 
one finds a contribution oc w^, 

ImE(*,.) = ^ ^ (=1)^ ( (-^f In r^^) . (36) 

which is seen to be less divergent than the contribution from even modes. 

In the conclusions we will discuss the effects of mass enhancement on 
thermodynamics and the feed-back to the collective modes. 
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5. Resistivity 

The contribution of critical fluctuations to the resistivity involves the 
solution of a quantum kinetic equation. For a Galilean invariant system the 
solution is simple: the center of mass of the charged particles accelerates 
freely in the presence of an electric field giving rise to a vanishing resistiv- 
ity. Only if the momentum is relaxed sufficiently rapidly due to impurities 
or Umklapp processes on the Fermi surface in a lattice system, one can 
calculate the resistivity from the transport relaxation rate. There are two 
contributions to this rate, from scattering-out processes and from scattering- 
in processes. The former is characterized by a relaxation rate which is twice 
ImE(efc,/c). The latter is an integral operator expressing the effect of ver- 
tex corrections. It is well known that for an isotropic system the vertex 
corrections reduce the probability for a collision in which momentum q is 
transferred by a factor 1 — cos^ = {q/2kp)'^. Considering that the typi- 
cal momentum transfer in a collision process with a critical fluctuation is 
Aq ~ this reduction factor is on average (1 — cos 9) ~ (fci?^)^^. 

Using these simple arguments, we expect for the transport scattering 
rate of quasiparticles 



— = [kFq J dek21mE{ef,k)l^-^J (37) 



and for the resistivity 



m* 1 , , 

Ap=^ . (38) 



With J de/s ImS(efc, A;) ^ ~ — lniTi{kFT), we obtain from Eq. (|26() in three 
dimensions (within the assumptions specified above) 



(39) 



where Tp is the bare Fermi temperature, i.e. unrenormalized by critical 
fluctuations (the noncritical renormalization is, however, included). It is 
seen that Ap increases linearly with the correlation length ^, and has the 
usual Fermi liquid dependence. In d = 2 dimensions we obtain 



A,.^^i <£|);(«)^(^)%„,r?/T) ,40) 



which depends quadratically on ^ and has a logarithmic correction to the 
temperature dependence. 
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6. Conclusions 

Following Ginzburg and Landau, quantum phase transitions are most 
often described in terms of order parameter fluctuations, i.e. bosonic collec- 
tive degrees of freedom. Recently, there has been an intense debate whether 
such a description is in principle possible in metals where always both bosonic 
and fermionic modes coexist at a quantum critical point, for a recent review 
see Ref. El The situation is much simpler in the regime considered in this 
article: in the low-temperature limit of the symmetric phase, a Fermi liquid 
description should always be possible, allowing for a controlled calculation 
of low-energy properties. It is, however, not possible to calculate all physi- 
cal quantities (like the momentum dependence of the collective modes and 
the mass renormalization) rigorously within the Fermi liquid approach as 
they are influenced by high-energy modes. Nevertheless, with a few extra 
assumptions on the regularity of the high-energy theory, we were able to 
estimate these terms reliably. 

Most of the results obtained in this paper are consistent with what is 
expected from the purely bosonic theory. For example, the specific heat 
coefficient 7 of a theory with dynamical critical exponent z is expected to 
be of the fornP 7 ~ T z at the quantum critical point and accordingly 
7 ~ away from criticality for T — > (with logarithmic corrections for 
d = z). This is precisely what we obtained within the Fermi liquid approach, 
see Eq. ((23) and (jSU). 

Using the controlled Fermi liquid calculation one can check which as- 
sumptions underlying a bosonic theory are valid. For example, for the 'even' 
modes we find that all Fermi liquid corrections to the order parameter suscep- 
tibility in p2|) and (|1I-{I14|) cancel as Np/m* and m*VF are not renormalized. 
A similar result was obtained by Chubukov, Pepin and ReclJ^who found that 
there are no singular self-energy corrections to x at the z = 3 critical point. 
Similarly, there are no singular Fermi liquid corrections to the self-energy 
of quasiparticles (|26l27j) besides the trivial Z factor, again consistent with 
literaturJ^. The situation is quite different for the 'odd' modes where the 
m* /m corrections do not cancel in Eqs. (|13j) and (|14j) . This implies a break- 
down of the Hertz approach for these modes and changes, e.g., the lower 
critical dimension of the Pomeranchuk quantum critical point which is an 
interesting topic for future investigations^^. 
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